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We show that a rectangle triangle random tiling with a tenfold symmetric phase 
is solvable by Bethe Ansatz. After the twelvefold square triangle and the eightfold 
rectangle triangle random tiling, this is the third example of a rectangle triangle 
tiling which is solvable. A Bethe Ansatz solution provides in principle an accu- 
rate estimate of the entropy and phason elastic constants. In the twelvefold and 
eightfold cases even exact analytic expressions have been obtained from the Bethe 
Ansatz solution. 



1 Introduction 

A random tiling ensemble with a tenfold symmetric phase is defined by rect- 
angles and isosceles triangles of sides of length 1 and I — 2cos(37r/10) = 
\/2 -\- T It, where X — (\/5+ l)/2 is the golden mean. This random tiling was 
used by He tt am and by Nissen and BeeliEl to model a decagonal phase in 
FeNb, by Oxborrow and |Mihalkovica to model disorder in decagonal AlPdMn 
and by Roth and Henleyu to model the equilibrium structure resulting from a 
molecular dynamics simulation. The perfect quasiperiodic tiling corresponds 
to a maximum possible density of a decagonal disc packing. B The aim of this 
work is to calculate the entropy of the random tiling and its phason elastic 
constants. To perform i.his calculation we use a transfer matrix that generates 
the ensemble of tilingsB The logarithm of the largest eigenvalue of this matrix 
will give the free energy in the thermodynamic limit. 

2 Bethe Ansatz 

We show that the model is solvable in the sense that its transfer matrix can 
be diagonalizcd using Bethe Ansatz techniques. First we deform the tiling 
such that its vertices are a subscLof those of the square lattice. Similar to t 
twelvefold square-triangle tilingLrO and the octagonal rectangle-triangle tiling 
we have to decorate the deformed tiles with lines to make this deformation map 
bijective. An example of the deformation and decoration is given in Fig.^. 
From this example it is seen that the triangles form domain walls between 
patches of one type of rectangle. There are three types of domain wall; one 
type is running from bottom- left to top-right (right mover) , denoted by the gray 
lines in Fig.Qb, the other two are running form bottom-right to top-left (left 



Figure 1: a) Patch of the random tiling, b) Corresponding patch on the square lattice 



movers) and are denoted by the solid and dashed hnes in Fig.|I|b. These hnes 
scatter with each other, but they persist through the lattice. This means that 
the number of each type of domain wall is a conserved quantity under the action 
of the transfer matrix. We denote the numbers of gray, solid and dashed lines 
by m, ni and n2 respectively. A state a on a row of the lattice can be specified 
by the positions yi, . . . , ym of the right movers and by the positions xi, . . . ,Xn 
of all the left movers n = ni+n2 with the specification that the lines «i, . . . , 
at positions Xi^, . . . ,Xi^^ are dashed lines. Elements ipio') of the eigenvector 
thus can be written explicitly as ipiH, ■ ■ • i l^^i: • • ■ i ^n', J/i: • • ■ , Um)- We make 
the following (plane wave) Ansatz for the form of the eigenvector 

. . . jirijxi, . . . ,a;„;2/i, . . .,y„i) = 

71 7n 7l2 ic 

j2 E ^(r)B(A^) n n < n n ^(^^^ ^^r). (i) 
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The amplitudes A depend on thcjpcrmutations tt and p and on the configura- 
tion of the various domain walls Simlarly, the amplitudes B depend on the 
permutation fi and on the sequence of dashed and solid black lines. 

If all the domain walls are separated the action of the transfer matrix is 
just a shift of each line to the right or to the left. The eigenvalue corresponding 
to the vector (|l|) is therefore given by 

n m 

^ = 11' IK (2) 

i=l 3 = 1 

At places where different domain walls are close together, the action of the 
transfer matrix is not given by a mere shift of all domain walls. These more 



complex processes turn out to be consistent with the form (|l|) of the eigenvector 
and impose constraints on the amphtudes A and B and the z's, w's and m's of 
the eigenvector (|l|). These equations are the socalled Bethe Ansatz equations 
and are given by 

rn . 
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fc=l ^ ' 4=1 



n f n (^^^^- + -r^^"7') n-(^^)' (4) 
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where u is given by 

wG,*)-(i;,+zf-zr2)-\ (6) 

The Ansatz (|l|) and the Eqs. can be generahzed to include chemical 

potentials to control the tile densities in the ensemble. 



3 Entropy 

The entropy can now be calculated from the eigenvalue (|2|). Some finite size 
results for the entropy per lattice site and extrapolation are given in Table ra. 



Table 1: Numerical data for the entropy 



N 


On 


8 


0.19482 


13 


0.17779 


21 


0.17088 


34 


0.16853 


55 


0.16763 


89 


0.16730 


oo 


0.1671 (1) 



According to the random tiling hypothesis the entropy per area has its 
maximum at the point where the densities match those of the quasiperiodic 
tiling. For the tenfold tiling the expression for the entropy in terms of the 



phason strain near the quasicrystalline point is given by 



aa(E) = <Jo - -Ki {{El, + E2y f + {Ely - E2x)^) 

- {{El, - E2yf + {Ely + E2,f) - {El, + Ely) .(7) 

The last term in (|^) is present if the height component in the third commen- 
surate direction is 'rough'. This can be checked by calculating the finite size 
corrections to the entropy. For a critical system these go as 

vrcu, / 1 \ 

(8) 

where c is the central charge and Vg is a geometrical factor. The central 
charge measures the number of independent height components in these kinds 
of models. Indeed, we find c = 3. 

Although Eqs. (|^)-(|^) in principle can be solved numerically for very large 
values of N we found that the structure of the solution for the groundstate is 
numerically complicated, which is why we as yet cannot push the systemsize 
N larger than 89. This also has prevented us so far to give estimates for the 
elastic constants. 
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